Magnetic field dependence of the low-temperature specific heat of the electron-doped 

superconductor Pr1.g5Ceo.15 CuQ 4 
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We remeasured the magnetic field dependence of the low-temperature specific heat of the electron- 
doped superconductor Pri.gsCeo.isCuCU (Tc = 22±2K) under a different measurement procedure. 
Under field-cooling, the electronic specific heat follows C e i(H, T) — j(H)T from 4.5K down to 1.8K. 
In the field range Hci < H < 0.5Hc2, the Sommerfeld coefficient ^y(H) is well fit by a power-law 
j(H) ~ H 1 ^ 2 . This result suggests that the pairing symmetry is d- wave-like at all temperatures 
below 4.5K. Our new measurement shows no evidence for the linear field dependence of'y(H) found 
previously at T — 2K. 

PACS numbers: 74.25.Bt, 74.20.Rp, 74.25.Jb 



The pairing symmetry is an essential input to under- 
stand the mechanism of superconductivity and as well 
other novel properties of superconductors. It is generally 
accepted that the hole-doped (p-doped) high-Tc cuprates 
have a d-wave symmetry, based on most experimental 
results and theoretical studies^. However, the pairing 
symmetry of the corresponding electron-doped (n-doped) 
cuprates remains unclear. Phase-sensitive^, ARPES 5 ' 6 , 
penetration depthZ£ and Raman measurements^ suggest 
a d-wave symmetry. However, point contact tunneling 10 , 
penetration depthiiii^ and specific heal*^ revealed some 
s-wave features in the optimally-doped and overdoped 
samples, and led to the suggestion of a d-wave to s-wave 
symmetry transition upon doping or lowering the tem- 
perature below AK. Transport and ARPES measure- 
ments have shown that the Fermi surface of the n-doped 
cuprates evolves from electron-like to hole-like upon dop- 
ing, with a two-carrier behavior around the optimal 
dopin g i . Recently, it was proposed that Fermi sur- 
face topology may influence the pairing symmetry which 
lead to a symmetry mixing and/or a symmetry transition 
upon doping9*i&. 

The magnetic field dependence of specific heat has 
been confirmed to be a valuable tool for detecting the 
bulk pairing symmetry, since the electronic specific heat 
C e i in the vortex state [Hex <C H <C H C2 ) scales 
with H and \/H for s-wave and d-wave superconduc- 
tors respectiveljiii. However, it has been pointed out 
that proper measurement conditions are essential to pro- 
duce these theoretical predictions— In this paper, we re- 
port our new specific heat measurements on high-quality 
crystals of the n-doped superconductor Pri ^Ceo.isCuO^ 
(PCCO). Under magnetic field cooling (FC), the field- 
dependence of the low temperature (T < A.5K) electronic 
specific heat has the power-law behavior C e i ~ H 5 
down to 1.8K. This d-wave- like feature is consistent with 
the recent phase sensitive tunneling 4 and ARPES 6 ex- 
periments. These results supplant our previous specific 
heat data measured on similar crystals 1 ^ and rule out the 
possibility of a temperature-dependent d-wave to s-wave 
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FIG. 1: (color online) Total specific heat of two 
Pr2- a; Ce a; Cu04 crystals (crystal I + II) at various magnetic 
fields under field cooling conditions. The field is applied par- 
allel to the c-axis of the crystals. Straight lines are the linear 
fits (see text). 



symmetry transition at optimal doping. 

The detailed sample preparation and experiment setup 
have been described in our previous paper—. The opti- 
mally doped Pr 2 _ 3; Ce a; Cu04 (a; = 0.15 ± 0.005 deter- 
mined from WDX measurements) crystals were grown 
by the self-flux growth technique. Here we selected two 
PCCO single crystals (crystal I and II), each with thick- 
ness of ~ 50/1771 and a mass of ~ 2mg. After oxygen 
reduction, both crystals show a superconducting tran- 
sition at T ~ 22K, with a full transition width of 4K 
from SQUID magnetization measurements. The specific 
heat measurements were conducted in a Quantum De- 
sign PPMS, with a home-made specific heat puck— The 
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magnetic field was applied parallel to the c-axis of the 
crystals. We measured separately the specific heat of 
crystal I, and of crystal I and II (to enhance signal to 
noise at low temperatures). As stated above, our new re- 
sults were obtained under field cooling conditions. That 
is, for each data point with the same magnetic field, the 
field was first applied at 30K, which is well above the su- 
perconducting transition temperature (Tc). The specific 
heat was then measured after cooling to a fixed tempera- 
ture between 1.8K and (see Fig.^l. This procedure 
is crucially different from the previous measurement—, as 
we discuss later. Results obtained from different crystals 
show similar field and temperature dependence. Down 
to our lowest measured temperature (1.8K), no Schot- 
tky upturn was observed, which simplifies the analysis of 
the specific heat. 

In the vortex state with H < Hc2i the low tempera- 
ture specific heat of superconductors can be written as 
C = G p h + C e i , with C p h and C e i the phonon and the elec- 
tronic contributions respectively. The phonon contribu- 
tion Cph is represented by (3T 3 well below the Debye tem- 
perature. The electronic contribution C e i has two parts, 
C e i = ctT 2 + j(H)T, where aT 2 is due to quasi-particle 
excitations of a d-wave superconductor [a = for BCS s- 
wave superconductors) and y(H)T is due to excitations 
around the vortex cores and their extended regions^. 
7 (if) is the Sommerfeld coefficient, which varies with the 
magnetic field. In our temperature range of measurement 
the aT 2 term was not detectable so the specific heat is 
represented by C = (3T 3 + -f(H)T. The specific heat of 
crystal I and II from 4.5K to 1.8K is shown in Fig. ^ 
For clarity, we plot C/T vs T 2 at only a few magnetic 
fields. The straight line fits are evident for temperatures 
down to 1.8K and fields up to 6T. We note that the 
slopes of these straight lines are about 40% larger than 
the actual value of (3 for PCCO— because of a contri- 
bution from a small unknown amount of the addenda 
Wakefield thermal compound which was not subtracted 
from the data. We verified that the specific heat of the 
thermal compound follows C cx T 3 in this temperature 
range. As a result, the Sommerfeld coefficient j(H) can 
be determined from the the vertical intercept in Fig. ^ 
With increasing magnetic fields, these straight lines are 
shifted upwards, indicating the increase of "f{H). 

Fig. El shows the value of -y(H) from QT to 6T. A 
nonzero y(H = 0) is clearly seen, followed by a non- 
linear increase of 7 (if) with field and a saturation be- 
havior near 6T. Finite 7(0) has been observed in many 
cuprate superconductors and its origin has been inter- 
preted differently 19,20,21 . In this paper, we focus on the 
field-dependent part. Volovik pointed out that for a d- 
wave superconductor the quasi-particle spectrum in the 
extended regions of the vortex cores has a Doppler shift 
due to the supercurrent induced by the magnetic field—. 
As a result, the electronic specific heat from these excita- 
tions follows a field dependence by j(H) — AyH, in the 
range T <C Tq and H < Hci- On the other hand, for s- 
wave (fully gapped) superconductors, the quasi-particles 
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FIG. 2: (color online) The field dependence of the Sommer- 
feld coefficient 'y(H) of Pr2- a: Ce a; Cu04 crystals (crystal I + 
II). The solid line is the power-law fit j(H) = 7' + AH 1/2 in 
the field range of 0.2T < H < 3T. Inset: Plot of j(H) vs 
(/j, H) 1/2 ; Solid line: the linear fit of 7(H) = 7' + AH 1/2 . 



are thought to be confined to the vortex cores. Then 
the specific heat is proportional to the magnetic field, 
7(H) = AH, because of the linear increase of vortices 
with field. However, in real s-wave systems the power-law 
7 (if) = AH a with l/2<a<lis frequently seen, which 
makes the interpretation more subtle— For example, for 
the s-wave superconductor Nb, a deviation from a linear 
field fit is clearly seen above a crossover magnetic field 
H* ~ 0.25ifc2— For p-doped cuprate superconductors, 
the -y/i? dependence has been reported for many crystals 
in support of d-wave pairing symmetryiSiSi. However, 
there are also a few exceptions that need to be clarified 19 . 

We attempt to fit our specific heat data to a d-wave 
power-law (a = 1/2). In the field range of 0.2T < H < 
3T our data can be well fit by ^(H) = 7' + A\/~H, with 
7' = 1. 64m J/raol K 2 and A = 1.92±0.1mJ/molK 2 T 1 / 2 , 
as shown in Fig. [21 In the inset we replotted ^{H) 
against iJ 1 / 2 , which shows a very good linear relation 
in the fitting range. By adopting the theoretical form 
A = ln{-^§^) 1/2 we obtain j n = i.lmJ/molK 2 , if we 
take the reasonable choices Hci ~ 6T and a « 0.7—. 
Our d-wave-like fit is supported by the following: 
1) The y/H fit is confined to the magnetic field range 
of 0.2T < H < 3T. The deviation of 7 (if) from sJ~H is 
evident in both the low field (H < OAT) and high field 
[H > 3T) limit, as seen in the inset of Fig.0 The field 
range for the VW behavior actually validates the argu- 
ment that the d-wave fit only applies to the vortex state 
well below iy^^iLiSiSS. Thermodynamic studies show 
that the H C i is close to 0.1T for Pri. 85 Ceo.i5Cu0 4 25 . 
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FIG. 3: (color online) The field-dependence of the specific 
heat of P^-^CexCuO^ crystals (crystal I + II) at tempera- 
ture T = l.&K. The solid line is given by j(H) = 7' + AH 1/2 . 

The levelling off of the electronic specific heat at around 
{.iqH = 6T suggests that the Hci is about 6T in the 
current system— 

2) Our fitting parameter 7' is only slightly smaller than 
7(0). It is possible that 7(0) is enhanced by a small 
trapped magnetic field, which is usually observed in our 
measurement system. 

3) The high-field Sommerfeld coefficient 7 (if > Hex) 
is expected theoretically to be 7' + 7^— For our theoret- 
ical fit (previous paragraph), 7' + j n ~ 5.7 / LmJ/molK 2 , 
which is quite close to the measured value of j(H > 6T) 
shown in Fig. [21 

In order to compare with the results in our previous 
paper—, we also performed separate measurements and 
analysis at a few fields (n$H < LIT) at our lowest tem- 
perature T = l.&Bffl. In Fig. AC{H)/T = (C(ff) - 
C(H = 0))/T is plotted against magnetic field, by sub- 
tracting the zero-field specific heat. AC(H)/T is equiva- 
lent to A'j(H) — 7(ff) — 7(0), since the phonon contribu- 
tion is independent of magnetic field. The nonlinearity 
between A'y(H) and field is clearly seen. A power-law fit 
A 7 (fi) = AVH with A = 2.08 ± 0.1m J /mol^T 1 ' 2 , as 



shown by the solid line, is reasonably good for H > 0.2T. 

The power law dependence we find between the specific 
heat and the magnetic field at low temperatures is clearly 
different from our previous paper, where the power law 
exponent a was found to increase from 1/2 to 1 as the 
temperature decreased from 4.5A to 2A— . The principal 
difference between these two works is the measurement 
procedure. In the previous measurement—, the sample 
was cooled from above T c under a constant field ~ 9T, 
but then the data were taken by sweeping the field down 
at a constant temperature. In this work, we cooled the 
system under a fixed magnetic field from above Tc to the 
measurement temperature. The data in Fig.^&nd Fig.0 
were obtained this way by starting from above Tc for 
each field value. We found that the measurement proce- 
dure in our previous paper may have caused the variation 
of a for two reasons. First, the specific heat measured by 
the field sweep procedure is larger at low fields than our 
new data. This suggests that there is flux trapping in 
the sample during the field sweep and consequently the 
actual field value in the sample is larger than indicated 
by the PPMS. Second, below T = 4.5K, the temperature 
of our heat capacity stage can vary during the measure- 
ment period by up to 2% at 1.8 A (decreasing to 0% at 
4. 5 A"). This could lead to an artificial increase of C e i 
which varies with field. In our new work, the scheme of 
data extrapolation to T — > (Fig. ^| or data interpola- 
tion to T — > 1.8A (Fig. |3 26 eliminates these problems. 
We think that our new data measurement procedure is 
more appropriate for comparing C e /(ff) with theory. 

In summary, we remeasured the specific heat of opti- 
mal electron-doped PCCO crystals. Our new data show 
that the Sommerfeld coefficient, 7 (if), is fit nicely by 
VW in the magnetic field range Hci < H < 0.5ffc2 
and the temperature range 1.8A < T < 4.5A. This re- 
sult rules out a d-wave to s-wave symmetry transition 
in this temperature range for optimal-doped materials. 
Our result is completely consistent with the recent phase 
sensitive measurements in n-doped cupratesi Khodel et 
al. proposed a possible d-wave to s- or p-wave symmetry 
transition upon doping—, due to the change of the Fermi 
surface topology. Further work with more over-doped 
crystals will be necessary to address this issue. 
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